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Abstract 



Resonance processes are common phenomena in multiscale (slow-fast) systems. 

In the present paper we consider capture into resonance and scattering on reso- 

K* , nance in 3-D volume-preserving slow- fast systems. We propose a general theory 

of those processes and apply it to a class of viscous Taylor-Couette flows between 



(^ ■ two counter-rotating cylinders. We describe the phenomena during a single passage 

'nI" , through resonance and show that multiple passages lead to the chaotic advection 

^^ I and mixing. We calculate the width of the mixing domain and estimate a charac- 

f^ ' teristic time of mixing. We show that the resulting mixing can be described using 

a diffusion equation with a diffusion coefficient depending on the averaged effect of 

the passages through resonances. 



1 Introduction. 

Appearing in a variety of settings, multiscale systems present quite a challenge for direct 
numerical simulations (see e.g. [10]). Thus, a development of analytical tools is rather 
important. For a wide class of low-dimensional multiscale systems (that are also called 
slow-fast systems) where some variables change with a characteristic rate of order 1, while 
the other variables change with a rate of order of £ -C 1 the averaging method can be 
used for approximate description of the dynamics (see e.g. [1]). A standard condition for 
applicability of the averaging method is that the frequency of motion of the fast variables 
is indeed of order 1 everywhere in the space. The violation of the above condition near 
resonances of the fast system may give rise to the following single resonance phenomena: 



first, separatrix crossings or, second, scattering on resonance and capture into resonance. 
The hydrodynamic flows involving the separatrix crossings were introduced in [2, 16] and 
later discussed from the theoretical, numerical and experimental perspectives (see e.g. 
[20, 21, 8]). The theory of separatrix crossings for such flows was developed in [19, 14]. 

The flows where the chaotic advection is caused by scattering on resonance and capture 
into resonance were introduced in [4, 5]. The later works include [15]. However, the 
corresponding theory, that gives a quantitative description of either single crossings or 
the long-time dynamics was not developed. 

There are two major objectives of the present paper. First, we propose a general theory 
of scattering on resonance and capture into resonance in volume-preserving systems. In 
our study we rely heavily on the close analogy between volume-preserving and Hamiltonian 
systems, for which the theory of scattering on resonance and capture into resonance was 
developed in [12]. Second, we would like to demonstrate that the resonance phenomena 
can be used in real settings to create flows with coexisting mixing and non-mixing (KAM) 
domains with the respective size and mixing rates easily and independently manageable. 
In the process, we want to emphasize that the presence of resonances is not a sufficient 
condition for chaotic advection and to illustrate some of the necessary ingredients for the 
resonance phenomena to create the chaotic advection. 

The structure of the paper is as follows. In Sections 2 to 4 we discuss a general theory of 
passages through resonances in near-integrable volume-preserving 3-D systems. In Sect. 2 
we introduce basic equations and study the unperturbed and the averaged systems. In 
Sect. 3 we consider the structure of a resonance and discuss the breakdown of the method 
of averaging. Section 4 contains a detailed description of the general properties of the 
streamlines' dynamics in the vicinity of a resonance. Scattering on resonance and capture 
into resonance are discussed in Subsects. 4.1 and 4.2, respectively. We show that both 
capture and scattering destroy the adiabatic invariance and we obtain expressions for 
jumps of adiabatic invariants during a single passage. As the magnitude of the jumps is 
very sensitive to initial conditions, in the case of multiple passages through resonance these 
phenomena can be treated as random processes and we derive their statistical properties. 
In Sect. 5 we apply the technique developed in the previous sections to a particular 
problem of mixing between two counter- rotating cylinders. In Sect. 6 we discuss the 
long time dynamics of streamlines and present results of numerical simulations. We show 
that multiple passages through resonance lead to the destruction of adiabatic invariance 
and chaotic advection and estimate characteristic time of mixing. And finally, we show 
that mixing can be described using a single diffusion-type PDE, in which the diffusion 
coefficient is provided not by the molecular diffusion, but by the averaged influence of the 
resonance processes. Thus, we call the resulting diffusion the adiabatic diffusion. 



2 Main equations. The unperturbed and the aver- 
aged systems. 

In the first part of the paper (Sects. 2-4) we present a general theory of passages through 
resonances in near-integrable volume-preserving 3-D systems. In other words, we study 
the fiows that are a superposition of an integrable base fiow and a weak perturbation. 
While the inclusion of the time- dependence adds to the variety of possible fiow configu- 
rations (see e.g. [4, 15]) and extends the range of admissible resonances, we will see that 
our (time-independent) setting is sufficient to observe most of the related phenomena and 
is better suited for analytical description. 

In the present paper we consider the base fiows belonging to a class of integrable 3-D 
fiows that possesses two invariants. Introducing these invariants as new variables we can 
write the evolution equations in the following form: 

A = 0, 

/2 = 0, 

In such fiows almost all the streamlines are closed curves (except for those that connect 
the stationary points). Each curve is defined by the values of Ji and I2 and we denote 
it as Tj^j,^. Such systems belong to the two actions and one angle class of dynamical 
systems, where a single changing variable being a coordinate along a streamline. 

A generic small perturbation causes the invariants of the base fiow to change and a 
result fiow has the following form: 

Ji = evi{Iij2,0), 

h = ev2ilul2,0), (1) 

e = u;ih,l2)+egih,l2,0). 

System (1) is volume-preserving if the right hand side is divergence-free. For < £ <^ 1 
fiow (1) possesses two time scales. The variable 6 is fast and the variables Ji and I2 are 
slow with a ratio of characteristic velocities being of order of e. System (1) is an example of 
so-called slow-fast systems. Therefore, in the first approximation we can average velocity 
field (1) over the fast variable 6. It is clear that this approximation is valid everywhere 
except for a small part of the space where ^ ^ 0. The averaging over 6 reduces (1) to 

/i = e^ f vi{h,i2,e) de, 
12 = e^ f V2{ii,i2,e) de. 

27r Jo 



For the sake of brevity, in (2) we used the same notation for Ji and I2 as in (1). System (2) 
is 2-D volume-preserving and, therefore, Hamiltonian and integrable. The Hamihonian 
$(/i, I2) is the flux of the perturbation vector v = (f 1, f2, g) across a surface S spanning 
the streamline of the unperturbed system T/^ /j. The quantity $ is an integral of the 
averaged system (2) and it is an adiabatic invariant of the exact system (1) (see e.g. 
[19, 14]). If 9 does not vanish anywhere on a streamline, $(/i,/2) would be conserved 
with the accuracy of order e over times of order 1/e (see c). Note, that the conservation 
of $ is a direct consequence of the validity of averaging. 

3 Structure of the resonance and the breakdown of 
the method of averaging. 

It is widely known in the theory of dynamical systems that passages through the regions 
where the frequencies of a given system satisfy the rational condition (in terms of (1), 
riiEVi + n2ev2 + n3{uj + eg) = 0, where n^'s are integer numbers) lead to chaotic dynamics. 
As the ratio of characteristic frequencies in 1) is of order of e, for the resonance to be of 
a low order we must have uj = 0. This is the only important resonance in the systems 
similar to 1). In the process of motion some streamlines pass through the region where 
^ ~ 0, that we call a resonance zone. In the next two sections we develop a general theory 
of resonance crossings in volume-preserving systems analogous to the one first developed 
in [11] and then adjusted in [12] for Hamiltonian systems. The equation uj[Ii,l2) = 
defines a 2-D surface in the original 3-D space (a curve on the (Ji, I2) plane), that we call 
the resonant surface, or the resonance, and denote by R. In the vicinity of a resonance 
the variable 6 is not fast compared with Ji and 12- Hence, we cannot expect averaged 
system (2) to approximate the exact system adequately. As a result, the value of the 
integral of the averaged system, $, may change in the process of a passing through the 
vicinity of the resonance. 

Qualitatively a single passage through the vicinity of a resonance can be described 
as follows. A characteristic streamline of the exact system approaches the resonant zone 
with the value of $ oscillating with a small amplitude, ~ e, near some value <l>~. When 
in the process of the motion it arrives to the resonant zone, it is either captured into the 
resonance, or crosses the resonant zone without being captured. (Actually, there is also 
some intermediate regime of motion in the resonant zone, but it occurs for a very small 
measure of initial conditions; we will not discuss it.) Phenomeno logically, the difference 
between the two regimes of motion in the resonant zone can be described as follows. In 
the case of capture, upon the arrival to the resonant zone the phase 6 switches its behavior 
from rotation {6 changes monotonically, a streamline makes full circles around the axis) 
to oscillation {6 changes between two values, 9 changes the sign twice during each period). 
In the case of crossing the resonant zone without capture, 6 changes the sign only once 
in the resonant zone. As we show below, the two phenomena have qualitatively different 



impact on the streamlines' behavior. After the passage through the resonant zone (and 
far from the resonance) the value of $ oscillates near some other value, <l>^, again with a 
small amplitude ~ e. 

In the rest of the present section we provide a qualitative description of these phe- 
nomena (without any proofs) and return to a detailed quantitative description in the next 
section. 

In the case of capture into resonance, upon the arrival to the resonant zone a phase 
point drifts for a long time (of order ~ 1/e) along the resonant surface. As a result, $ 
changes by a value of order 1 . Among all the streamlines that arrive to the resonant zone 
during a given time interval of order ~ 1/e only a small part, of order ~ ^/e, is captured. 
Initial conditions for streamlines that are or are not captured are mixed. Therefore, it is 
reasonable to consider capture as a probabilistic phenomenon. 

The streamlines that cross the resonant zone without being captured typically pass 
through this zone in time of order ~ 1/v^ (see [11, 12] for more accurate estimates 
for Hamiltonian systems). The major resonant phenomenon for such streamlines is the 
scattering on resonance. In this case the difference, A$ = $+ — <l>~, that is typically of 
order ~ ^/e, is referred to as the amplitude of the scattering or, as $ is the adiabatic 
invariant of the averaged system, as a jump of the adiabatic invariant. This value is very 
sensitive to small changes of the initial conditions: a change of the initial conditions far 
from a resonance by a quantity of order ~ e leads to a big (of order ~ 1 compared with the 
typical values) relative change in the amplitude of the scattering. Hence, the scattering 
can be considered as a random process. 

As the passages through the vicinity of a resonance lead to the destruction of adiabatic 
invariance of $, the structure of streamlines becomes chaotic. Note, that the expectation 
of the change of $ (the probability of a phenomenon times the size of a change of $) 
is of the same order for scattering and capture. Thus both processes contribute to the 
stochastization equally. There is a crucial difference between these resonant phenomena 
and certain other routes to chaos, like separatrix splitting. Although the total size of the 
resonant zone is small with the magnitude of the perturbation and it is localized in the 
space, the effect is large and global: the chaotic domain is of the order 1 (of the order of 
the total size of the system) regardless of the magnitude of perturbation. 

4 The description of motion in the vicinity of the 
resonance. 

Let us derive equations of motion in the vicinity of the resonant surface, that is called a 
resonant zone. Apply change of variables 



such that 

da du du do- n j \ io\ 

where ui is the angular velocity defined above, a is an additional variable (chosen in such a 
way that (3) is satisfied) and /u(/i, I2) is a coefficient in the expression for the infinitesimal 
volume: dV = fi{li,l2)dlidl2d6. Actually, (3) can be relaxed to be required on R only. 
In terms of the new variables equations of motion (1) can be written as 

uj = efi{iu,a,e), 

a = ef2{uj,a,e), (4) 

9 = u + eg {u, a, 6) . 

and the divergence-free condition has the most simple form: 

duj da 89 
In the resonance zone we have: 

9 = ehiiu,a,9) + e^9 + Oie^) 



e h (^, ^, e) L=o + ^(^J + ^^]+0{e 
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It follows from the last line in the above equation that the width of the resonant zone is 



given by 



< consti/e. Rescaling time, 

r = A/e t, 



and denoting the derivatives with respect to r by a prime, in the leading order (4) can be 
written as 

cr' = v^/2,0, (5) 

e" = h,o-V-e^0'- 

In (5), we used a notation /i^o = /i (0, a, 9) and /2,o = /2 (0, (x, 9). We can expand fi^ as 

/i,o = a{a) + b{a,9), 




(a) (b) 

Figure 1: Schematic phase portraits on the {9,9') plane. 

where a{a) is the average value of /i^o over a period: 
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/i,o d9, 



and b{a, 9) is the remaining, zero-average part. Introduce a potential energy for /i_o: 

dV {a, 9) 



/i,o 
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The function V contains a periodic and a linear in 9 parts: 

V = a{(r)9 + I b{(r, 9) d9 = a{(r)9 + h{a, i 



Define the resonance energy: 



H, 



V\2 



+ v 



(6) 



and denote by h the value of Hr. Note, that at the crossing 9' = 0. 

The relation between a and max \b{9)\ defines the properties of the phase portrait on 
{9, 9') phase plane. If 

max \h{9)\ > \a\ , 

the phase portrait looks like the one shown in Fig. la. For the sake of simplicity, we 
assumed that there is only one oscillatory domain. The separatrix E is denoted by the 
dashed hne. In the opposite case, the phase portrait looks like the one shown in Fig. lb. 

4.1 Scattering on resonance. 

The value of $ changes in the process of motion near and inside the resonant zone. The 
cumulative change in $ during a single passage can be considered as a jump associated 



with a resonance crossing. One can see that the theory of scattering in volume-preserving 
systems is the same as in Hamiltonian systems developed in [12]. We refer the reader to 
that paper (see also [17]) for the details of derivation and the error estimates. 

Consider a streamline that crosses the resonant surface at some time moment t*. Along 
this streamline a phase point passes through the resonant zone fast enough, in time of 
order I/a/^. While the phase point moves inside the resonant zone changes of the slow 
variables are small, of order y^. Therefore, for approximate description of dynamics inside 
the resonant zone we can fix the value of slow variables at their values at the crossing. 

Let ti and t2 be two time moments, ti < t* < t2 such that |ti,2 — t^\ ~ 1/e and t* is 
the only moment of crossing of resonance on the time interval (ti,t2)- Then the jump of 
$ on the resonance is 

A$ = / $ dt. 

Jti 

In the first approximation, we can replace in the integral the limits t2,i with ±00 (far 
from the resonance $ does not change much). We get 

A$= / $rft = 2v^/ de = -2s^/e -M (7) 

y^oo J-soo e 9' y_,^ v^2 {K - V) e 

where s = sign(a) and 6^ and h^ are the values of 6 and h, respectively, when the exact 
trajectory hits R. 9^, is the resonance phase, that is equal to the shifted value of 9: 9^, = 
9^: — 2nK, where the value of K is chosen in such a way, that < —9^, — (b/a) sin ^* < 27r. 
If the quantity $ can be calculated explicitly as a function of Ji and I2, we can directly 
substitute $ = $(/i, I2, 9) into (7) and obtain the expression for A<I>. This is the path we 
take in the subsequent part of the paper, when we consider a particular example of the 
resonance-induced chaotic advection in Stokes fiows. However, in many problems $ can 
only be obtained in a linearized form in the vicinity of a resonance. In the latter case we 
use the change of variables (3) to write (7) as 

A* . _2.Vi r -=1= (^A„ . ^A,o) ^0. ,8) 

J-soo ^2 {K - V) \ du da J 

Define the quantity ^ as 

Zn \a\ 
where the curly brackets denote the fractional part. In terms of ^, (8) can be written as 

A$ = ~2sy^^= / , ^ = ^ V i.o + ^ V 2.0 d9. (9) 



j2|27r^ + ^ + 6(^)/a| 



The quantity ^ characterizes a scattering. For every set of initial conditions, the values 
of C, and A$ can be calculated exactly. However, a small change of order e in the initial 
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conditions produces in general a large change of order 1 in ^. Hence for small e it is the 
best to treat .^ as a random variable, uniformly distributed on (0, 1) (see [12]). 

Statistical properties of the scatterings depend on the shape of the phase portrait on 
the (^, 9') plane. If the phase portrait looks like the one shown in Fig. la, in a generic 
case, scatterings on resonance cause $ diffuse and drift. The rate of drift is determined 
by the ensemble average of A$ (the derivation is similar to [12], see also [13]): 

(A<|.) = y" A$(Orfe = -v^ ^^i^i^(0' ^) ~ v^' (10) 

where \E'(0,o") is the flux of /2,o through the area (that we denote by Sr) inside the 
separatrix S on the {9, 9') phase plane in Fig. la: 

^{0,a)= f f f2,od9d9'. (11) 

Note, that in Hamiltonian systems (and in many volume-preserving systems, in particular 
in the system discussed below) f2fl does not depend on 9 and, therefore, /2,o can be taken 
out of the integral and \E' becomes proportional to Sr. 

If the phase portrait looks like the one shown in Fig. lb, the mean change of $ is zero: 
as there is no separatrix, \1'(0, a) = 0. 

The rate of diffusion depends on the dispersion of A$: 



D= [ {A^O-imfd^r^e. 
Jo 



'0 

We return a more detailed description of the long-time dynamics in Sect. 6. 

4.2 Capture into resonance. 

The other phenomenon that affects the behavior of streamlines at a resonance crossing is 
capture into resonance. 

Capture is possible only if the phase portrait in the {9, 9') plane looks like the one 
shown in Fig. la and the phase space (^, 9') has a well-defined separatrix. 

Capture into resonance can be described as follows. The flux \E'(0,o') changes as a 
phase point moves along a streamline. Suppose, \E'(0,o") increases. Then if a streamline 
comes very close to the hyperbolic fixed point, it may cross the separatrix and, as a result, 
be caught in the oscillatory domain within the separatrix loop. In this case, the streamline 
starts moving on the resonant surface. 

Captured dynamics consists of fast rotation on the (^, 9') phase plane and slow evolu- 
tion of the parameters of the orbit, a and h. The slow equations are 

a' = v^ (/2,o), (12) 



where (■) denotes the averaging over the fast, {0,0'), motion. The slow dynamics is 
governed by a new conservation law: the flux of /2,o through the area inside a closed 
curve, that we denote by rij(/i, o"), on the {0,9') phase plane is constant: 



m{h, (y)= /2,o dO dO' = const. 

J JrR{h,a) 

Recall, that the curve rj:j(0, a) corresponds to the separatrix and we come back to defini- 
tion (11). The captured motion is integrable and \l'(/i, a) is the Hamiltonian. 

While the streamline moves, the flux through the separatrix loop, that is a function of 
slow variables, changes and, if it returns to the initial value, the streamline is released from 
the resonance. If the flux keeps increasing, the once captured streamline stays captured 
forever (in real settings, until it reaches the limits of the system) . 

As it was discussed in Sect. 3, capture can be considered as a probabilistic phenomenon: 
initial conditions for streamlines that are or are not captured are mixed. Consider a 
point M far from the resonance such that streamline passing through M intersects the 
resonance. Let V^ be a sphere of radius 6 centered at M and V^'"^ be the part of V^ formed 
by initial conditions of trajectories with a capture into the resonance (see [11, 12, 9]). We 
define the probability of capture for the streamlines starting at M as 



P{M) = v^ lim lim '^ /\ . 

^ ^ 5-.0 e->0 vol V^ 



Following [12], we have: 



P(M) = V~e^^^P^^I^^-V~e, ii{dSR/da)^>0, (13) 

I I * 

where the subscript '*' indicates that the corresponding quantity must be evaluated at 
the resonance. For {d'^{0,a)/da).^ < 0, P{M) = 0. 
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5 An example: Taylor- Couette Stokes flow. 

In the rest of the paper we discuss the chaotic advection and mixing in a particular 
example of volume-preserving system. Besides illustrating a general theory developed in 
the previous sections, we would like to demonstrate that the resonance phenomena can 
be used in real settings to create flows with coexisting mixing and non-mixing (KAM) 
domains with the respective size and mixing rates easily and independently manageable. 
As a model base problem we chose a relatively simple system - a Taylor-Couette flow 
between two infinite counter-rotating circular coaxial cylinders. Following [5], we restrict 
our consideration to Stokes (small Re) flows. However, we would like to stress that the 
whole phenomena is structurally stable, which means that it survives most of the changes 
in the geometry of the flow settings and peculiarities of the flow perturbations. 

In the system under consideration, the unperturbed flow is given by (a solution to a 
2-D Stokes equation): 



r = 0, 

z = 0, (14) 

e = uj{r,z), 



with 



U{r, Z) = UJi— ( 1 - ^ ( 1 - T]— ] ] + iOi— ( 1 - T] — 



Ti \ 1 — r]^ \ ^i J J ^ \ oOi J 1 — r]"^ 

where 

V = n/ro- 

In the above equations, Tj, Tq are the radii of the inner and outer cylinders, respectively. 
In what follows, we assume that the inner cylinder rotates with a constant angular speed 
UJi, while the angular speed of the outer cylinder changes periodically with z around the 
average value of — cjj: uJo = cjj(— 1 + 6sin{Xz)), where A and 6 are the wavenumber and 
amplitude of oscillations, respectively. Qualitatively same results (differing in the exact 
expression for the adiabatic invariant and, consequently, in specifics of the resonance 
equations) can be obtained if, instead of the varying Uo, the outer cylinder has a varying 
radius. Similarly, we could consider a flow inside a torus-shaped container instead of the 
annulus. 

To introduce dimensionless variables, we rescale the time by uJi and the distances by 

r 
t = tiJi, p= —. 
r-i 

From now on we assume that all the variables are in the dimensionless form and the dot 

denotes the derivative with respect to t. For the numerical simulations presented below 

the following set of parameters was used: 

n = 0.2, To = 1.0, r] = 0.2, S = 0.4, A = 27r. 
11 



In what follows, we assume periodic boundary conditions in (normalized) z with a period 
1. 

The value of p changes between p = 1 (at the inner cylinder) and p = l/rj (at the 
outer cylinder). For uj{p,z) we have: 

T7 1 1 7] / 1 \ 

^(P, z) = -P-, + --, + -, 2^ sm{\z) P - - • (15) 

1—1] pi — 7] 1—1]^ \ pj 

A perturbation consists of two parts. First, the inner cylinder is propagating upward (in 
the axial direction) with a constant (slow) speed z = e. Second, the shape of the outer 
cylinder deviates slightly from being exactly circular. The full 3-D flow is given by 

p = e/t (p — 1) cos^, 

z = £ (1 + Inp/lnr]) , (16) 

6 = uj(p,z) eK,(2p — l)sm6. 

P 

In (16), < e <^ 1 is a small parameter, while /t ~ 1 defines a characteristic ratio of 
the two perturbations. By changing the value of k, we can vary the properties of the 
resonance phenomena (see Subsect. 5.3 below). One can see that the axial velocity, z, 
equals e at p = 1 and vanishes at p = l/rj. A limiting case k, = corresponds to the outer 
cylinder being circular and the axial symmetry of the flow. Note, that the reason for 
choosing both perturbations to be of the same order is that it is the interaction between 
the perturbations that results in the chaotic advection. 

5.1 The averaged system and the structure of the resonance. 

One can see that in (16) the variables p and z are slow and the variable 6 is fast. In the 
first approximation we can average (16) over 6 to get the following averaged equations of 
motion (cf. (2)): 

p = 0, 

(17) 
z = e {1 + lnp /Inrj) . 

The averaged trajectories (in the full 3-D, {p,z,6), space) wind the cylinders of con- 
stant radius with the direction of the rotation depending on the value sign of uj. System 
(17) is Hamiltonian with Hamiltonian 

6$ = er, (p+:r^{p Inp - (p - 1))) . (18) 

\ In rj J 
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The quantity $ is analogous to H in [5]. It is an integral of the averaged system and 
is an adiabatic invariant of the exact system: in the absence of resonances it would be 
conserved with the accuracy of order e over times of order 1/e (see [3]). 

The averaging is valid away from a resonance surface R, where tu = 0, that is given by 



Pr{z) - 



1] 1 + 1] — 6sm{Xz) 

Trajectories that do not intersect R remain regular. As in the averaged system the value 
of p does not change (see (17)), in the first approximation, the chaotic domain occupies 
the annulus between pmin and pmax, that are the minimum and maximum values oi pfi{z), 
respectively, and are given by 



1 1 + 1] + ri5 11+1] — ri5 

Pmin \ I ! | ^ anCl Pmax \ I ! ^ • 

1] 1 + 1] + y 1] 1 + 1] — 

For the above values of the parameters, we get that the mixing domain is located between 

Pmin = 2 and Pmax = V7 ^ 2.65. (19) 

Note, that pmin and pmax do not depend on the strength of perturbation, e. The division 
of the fiow domain on the chaotic and regular subdomains is shown in Fig. 2. The parts 
to the left and to the right of the corresponding vertical lines are filled with KAM tori. 
One can see that for 5 = 1, we get pmax = 1/v (i-e. rmax = ro), in other words, all the 
outside part of the annulus is chaotic. To mix the part adjacent to the inner cylinder, it 
is more effective to modulate the frequency of the inner, not of the outer, cylinder. In the 
opposite limit, as 5 -^ 0, the chaotic domain shrinks to the zero width around p = a/1/^- 

5.2 Resonance variables. 

We can define a new variable a: 

a = -i^ r Pf^'\ , dz, (20) 

21] J, l-^sin(A^) ^ ^ 

Note, that a is just an auxiliary variable. All the final results will be given in terms of 
the original variables: p and z. Condition (3) is satisfied on R only. As o" is a function 
of z only, we can simplify the diffusion equations below. In the right hand side of (5) we 
have: 

/i,o = a + bi cos 9, 

1] f 1\ 
" = -6\ cos{\z) [ p ) (1 + 1]^P /l]^ v) 



1 — 1]^ \ p 

1 



hi = -2k 

P + 1 
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Figure 2: Division of the flow domain. A chaotic domain is between the vertical lines. A 
regular (KAM) domain consists of two parts at the left and at the right. The wavy line 
in the middle is the resonance. 



and 



f2,o = ~p{p'-l) (1 + lnp/lnr/). 



The relation between a and 6i defines the properties of the phase portrait on {6, 6') phase 
plane. If 

\bi\>\a\, (21) 

the phase portrait looks like the one shown in Fig. la. If (21) does not hold, the phase 
portrait looks like the one shown in Fig. lb. Note, that a = at A2; = 7r/2 and Xz = 37r/2 
(where cos(A2;) = 0). Therefore, while in the center of the mixing region condition (21) 
may be either satisfied or not, there are always zones near the edges of the mixing regions, 
where (21) is satisfied. 



5.3 Scattering. 

In the problem under consideration, we have an explicit form of $ = $(p). Hence, we can 
use equation (7) to calculate the jumps A$ (or Ap). We have: 



A$ 



^ dt = -2s\fe TiK (p - 1) (1 + Inp /In i]) 



cos 9 



d9. 



(22) 



„ y/2 {K - V) 

where s = sign (a). The value of p must be taken at the moment of crossing, p = p^,, and 
the resonance potential, V, is 

V{a,e) = -ae-bisme. 

In terms of C,, we can write (22) as 

(p-1) r^* cos 5 



A$ 



-2sv^i 



SK r,- 



[1 + Inp /In 7]) 



^/2\27r^ + e + {bi/a)sme 



d9. 



14 



with 

^ = {:z\-A = ^\-e.-ih/a)sme,\. 

In \a\ In ' ' 

The quantity 9^ was defined in Sect. 4.1. Statistical properties of the scatterings depend 
on the shape of the phase portrait on the {6, 9') plane. If the phase portrait looks like the 
one shown in Fig. la, the ensemble average of A$ is (see Sect. 4): 

(A$) = -2s^e Kr.^A^ (1 + Inp /In v) Sr = s^e ^— ^ (1 + Inp /In r^) Sr, 
Znbi Zn 

where Sr is the area under the separatrix loop in Fig. la: 

Sn = 2 f ""' y/-2{V - V,)d9, 

^ "rain 

and y^ is the value of V at the hyperbolic fixed point in Fig. la: 

Vc = —a9c — h\ sin^c; cos^c = —a/hi. 

If the phase portrait looks like the one shown in Fig. lb, (A$) = 0: as there is no 
separatrix, Sr = Q. 

As $ is a function of p only (see (18)), we can use p instead of $. We have: 

r- P — 1 f^* COS 9 ,„ , , 

Ap = -2sJeK^-^ - M 23 

^ ' ^2|27r^ + ^ + (6i/a)sin^| 



Equation (23) was checked numerically for various values of parameters ^, k, and e. In 
Fig. 3 the plots of Ap (^) /^/e are presented for (a) /t = 2 (when (21) is satisfied) and (b) 
K = 0.2, (when (21) is not satisfied). The solid lines in Fig. 3 correspond to theoretical 
values of Ap{e) /y^ given by (23), and the asterisks show values obtained numerically 
from (14) for various values of ,^. When inequality (21) is satisfied, Ap{^) has a singularity. 
Therefore, there is a possibility (albeit quite small) of large changes in p in the process of 
scattering. However, as this singularity is logarithmic, (Ap) is finite: 

(Ap) = -sV^ ^^Sn. (24) 

Zn 

The plot |(Ap)| /^/e as a function of p for k, = 2 (the values of other parameters defined 
in (19)) is shown in Fig. 4. Note, that the value of Sji (and, therefore, of (Ap)) depends 
crucially on the k and 6. In particular, S"/? — * as /t — * 0. 

Note, that the values of p in the right hand side of (23) and (24) are, strictly speaking, 
those at the moment of crossing. But, as characteristic values of Ap at a single crossing 
are small (it is the accumulation of those changes that leads to the mixing), if we specify 
the initial conditions far from R, we can use those values of p in (23) and (24). 
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Figure 3: The plot of Ap/y/e as a function of ^; (a): k = 2 and (b): k = 0.2. Note the 
difference in scales. 
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Figure 4: The plot of |(Ap)| /y^ as a function of p for /t = 2 and the values of other 
parameters defined in (19). 

It follows from the definitions of s and a, that (Ap) > if cos(A2;) < and (Ap) < 
if cos(A2) > 0. As every single jump is small, two consecutive crossings occur at almost 
the same values of p and Sji. Therefore, it follows from (24), that the average change in 
p during one 2;-period is zero. 



5.4 Capture. 

A small percentage of particles do not leave the vicinity of resonance quickly, but are 
captured into resonance and travel near R for a long, of order of ~ 1/e, time. 

It was discussed in Sect. 4, that capture can be considered as a probabilistic process 
with a probability to be captured from a small ball of initial conditions given by (13): 



P{M) 



27r lal 



2-7? lal 



R 



dS_ 
dz 



dz 



(25) 
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Figure 5: Captured motion, (a) A projection of a streamline on the slow, (p, -z), plane; 
and (b) the time evolution of uj{p, z). 

A typical captured dynamics (for e = 10~^ and k = 2) is shown in Fig. 5. A projection 
on the slow, (p, z), plane and the time evolution of a;(p, z) are shown in Figs. 5a and 5b, 
respectively. A streamline comes from the bottom in Fig. 5a (from the left in Fig. 5b), 
is captured near z = 0.05 (t = 100), moves along the resonance, is released from the 
resonance near z = 0.45 (t = 1000), and then proceeds along an adiabatic path. To 
describe the captured dynamics we could use equations (12) with expressions for /2,o and 
V derived in the previous subsection. But as /2,o does not depend on 9, the evolution 
equation for a, 

Cr' = V^/2,o(o-), 

can be solved separately. As a depends on z only, we can return to the original variables 
to get 



v^(l + lnp/lnr/)^^o = v^ 1 + 



1 



l,ui-' + " 



ri6sm{\z] 



21nr7 \ri 1 + rj — 6 sin{Xz) 



(26) 



Note, that equation for the captured motion (26) can be solved once regardless of the 
initial conditions. Different initial conditions (i.e. different moments of capture) can be 
taken into account by shifting the origin of the slow (captured) time. The averaged (over 
fast {0,9') oscillations) evolution p{t) is given by the resonance condition uj{p,z{t)) = 0. 
It follows from the symmetry Xz ^ it — Xz that the release occurs at approximately the 
same value of p (with the accuracy of order of y^) at which the capture happened. It 
takes a phase point time (recall, that the original time differs from the slow time r by a 
factor l/y/e) 

T, = ^ / -dz 

V ^ Zcap 

to travel from capture to release. 
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Figure 6: Regular and chaotic domains. Almost straight vertical lines are regular stream- 
lines. A single streamline that starts at rj„ = 0.45 fills almost the entire chaotic domain. 

6 Long term dynamics and numerical simulations. 

It was shown before (see e.g. [17] and references therein) that in the Hamiltonian systems 
the accumulation of jumps of an adiabatic invariant at resonances leads to the chaotic 
advection and mixing. In the present section we study the dynamics over the long intervals 
of time that include many (of order of ~ 1/e) crossings. 

There are two quantities that describe the chaotic advection: the size of the chaotic 
domain and a characteristic rate of mixing inside the chaotic domain. 

It was shown in Sect. 5, that the chaotic domain is a cylinder between Vmin and Vmax- 
Outside the chaotic domain the majority of streamlines are regular. More precisely, they 
are regular if they are removed from the resonance by a distance that is small with e. A 
technique to estimate the excess width of the resonance domain was developed in [18] . A 
projection of three representative streamlines on the slow, (r, z), plane is shown in Fig. 6. 
Almost straight vertical lines are regular streamlines. A single streamline that starts at 
Tin = 0.45 fills almost the entire chaotic domain. 

Note, that the size and the location of the chaotic domain depends on the size of the 
annulus, rj, and the amplitude of the frequency oscillations, S, and is independent of the 
magnitude of perturbation, e. The size of the chaotic domain is of order 1 (i.e. it is on the 
scale of the whole system) regardless of how small e is. This property sets the resonance 
phenomena (and separatrix crossings, see e.g. [19]) aside from most of the other "routes 
to chaos." While the phenomena itself is local (we need to solve equations of motion only 
in the vicinity of a resonance surface), the effect is global and does not disappear as e 
becomes smaller and smaller. 

Recall that for £ = the dynamics is regular (see Sect. 5.1) and we may expect a 
smooth transition from partially chaotic to fully regular dynamics as e — i> 0. The only 
chance to achieve it is if as e — ;> it takes longer and longer time for streamlines to mix. 
In other words, the rate of mixing (diffusion) must go to zero, as e ^ 0. To estimate the 
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rate of diffusion we need to know tlie statistics of the changes of the adiabatic invariant. 
In a generic case, depending on the ratio between the coefficients in V (the resonance 
potential), the value of $ may drift or diffuse. However, in the particular case of the flow 
given by (14), it follows from the symmetry Xz ^ n — Xz, that (Ap) at two consecutive 
crossings cancel each other. Therefore, the statistics of the changes of p is diffusive 
regardless of the validity of (21). 

Quantitative properties of the diffusion of adiabatic invariant depend on whether con- 
secutive crossings are statistically dependent or independent. Statistical independence 
follows from the divergence of phases 6 along trajectories and for Hamiltonian systems 
was discussed in [12]. In volume-preserving systems it can be illustrated following the sim- 
ilar lines (see also [6]). Assuming the statistical independence of consecutive crossings, we 
can estimate the rate of mixing. The evolution of p can be considered as a random walk 
with a characteristic step of order of ~ ^/e. Hence, after A^ crossings a value of p changes 
by a quantity of order of ~ viV x ^/e. The mixing can be called complete if a trajectory 
spreads all over a chaotic domain. The difference between the values of p that bound the 
chaotic domain in our problem is of order 1. Therefore, it takes approximately N ~ 1/e 
resonance crossings to complete the diffusion. As a typical time between successive cross- 
ings is of order of ~ 1/e, a characteristic time of mixing, Tm, is of order of ~ e^'^. One 
can see that, like we suggested above, a characteristic time of mixing goes to infinity as 
e ^ and the rate of mixing, defined as 1/Tm, vanishes as e ^ 0. The similar estimate 
was reported in [7]. 

Another important limiting case corresponds to k ^ 0. For k, = the outer cylinder 
is circular and we regain the axial symmetry. In this case the resonances do not lead to 
chaotic advection. Indeed, it follows from (23), that a characteristic size of the jumps 
Ap ^ as K ^ 0. 

To check to the validity of theory developed in the previous sections, we performed a 
set of numerical simulations. We used the set of parameters specified in the beginning of 
Sect. 6 and 

e = 10"^; K = 2. 

First, we checked the validity of formula for a jump in adiabatic invariant at a single 
crossing (Eq. (23)). Those results were presented in Sect. 6.3. Second, we studied the 
diffusion of adiabatic invariant and large scale mixing. For that purpose we took 1000 
initial conditions that are uniformly distributed in a box of pin x Zin x Oin = [2.249, 2.251] x 
[—0.01,0.01] X [—0.01,0.01] and considered the Poincare sections located aX z = A^ -1-0.25 
and z = N + 0.75, where A^ is a set of integer numbers. Every trajectory crosses the 
resonance once between the consecutive sections. However, to simplify the following, it is 
convenient to consider a change in p after two successive crossings, Ap^ = Apo.25-^0.75 + 
Apo.75_»o.25, as Aps has zero mean. 

Let \E'(p, A^) denote the number of trajectories that, after A^ double crossings, have 
the value of p between p — 0.005 and p + 0.005. The spreading of ^(p, A^), obtained by 
integrating (14) with the initial conditions specified above, is shown in Fig. 7. 
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Figure 7: The histogram of \l/(p, A^) after different numbers of the double resonance 
crossings. The sohd hne and the stars correspond to A^ = 10 and A^ = 400, respectively. 
The other curves are between those values of N. 

Many systems with random walks are described using diffusion equations, that for a 
system under consideration can be written as 



d^ 



D 



dt 



d-p r ^^^ -w 



(27) 



where "^n is a probability distribution function. To relate Eq. (27) to the evolution of 
\l/(p, A^), the time t should be measured in the units of the period in the z-direction 
(in other words, there is one double resonance crossing per unit time). The diffusion 
coefficient -D(p), that we call the coefficient of the adiabatic diffusion, is given by the 
dispersion of Ap^: 



D = D{p; 6, K, A) 
In the leading order we have 



{Ap{0-{Ap)fd^. 



D{p] e, K, A) = eD{p] k, A). 

The profile of the normalized dispersion, D, as a function of p for the values of the param- 
eters specified above is presented in Fig. 8. The second moment of \l/(p, A^), is presented 
in Fig. 9. The constant slope confirms the diffusion assumption and the magnitude of 
the slope is in the good agreement with the prediction based on the simplified diffusion 
equation with D{p) = D{p = 2.25) ^ 5 x 10"^. 

The above considerations describe the diffusion across the level of constant adiabatic 
invariant (in radial direction). Mixing in the axial direction is due to the radial gradient 
in the axial velocity, while mixing in the azimuthal direction is due to the gradient in 
uj{p, z). Both processes are much faster then the radial diffusion and by the time of order 
of 1/e^ the distribution in those directions is quite uniform. 
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Figure 8: The normalized dispersion, D, as a function of p. 
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Figure 9: The square of the standard deviation, ((p — po) ); averaged over 1000 trajecto- 
ries as a function of the number of the double resonance crossings, A^. 

Conclusions and Acknowledgements 

In the present paper we considered capture into resonance and scattering on resonance 
in 3-D volume-preserving slow-fast systems. In the first sections we proposed a general 
theory of those processes. Then we applied it to a class of viscous Taylor- Couette flows 
between two counter-rotating cylinders. We described the phenomena during a single 
passage through resonance and showed that multiple passages resulted in the chaotic 
advection and mixing. We calculated the width of the mixing domain and estimated a 
characteristic time of mixing. We showed that in the systems with resonances mixing 
can be described using a diffusion equation with a diffusion coefficient depending on the 
averaged effect of the passages through resonances. Finally, we would like to stress that 
the resulting diffusion coefficient is not related to the molecular diffusion. 
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